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Abstract: Optical absorption is omnipresent and very often distributed 

non-uniformly in space. We present a numerical study on the effects of 
inhomogeneous absorption on transmission eigenchannels of light in highly 
scattering media. In the weak absorption regime, the spatial profile of a 
transmission channel remains very similar to that without absorption, and 
the effect of inhomogeneous absorption can be stronger or weaker than 
homogeneous absorption depending on the spatial overlap of the localized 
absorbing region with the field intensity maximum of the channel. In 
the strong absorption regime, the high transmission channels redirect the 
energy flows to circumvent the absorbing regions to minimize loss. The 
attenuation of high transmission channels by inhomogeneous absorption is 
lower than that by homogeneous absorption, regardless of the location of 
the absorbing region. The statistical distribution of transmission eigenvalues 
in the former becomes broader than that in the latter, due to a longer tail 
at high transmission. Since the maximum transmission channel is the most 
efficient in bypassing the absorbing region, the ratio of its transmittance to 
the average transmittance increases with absorption, eventually exceeds the 
ratio without absorption. 
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1. Introduction 

The ability to control light propagation in turbid media is of great importance to many fields, 
ranging from medical imaging, laser surgery to photovoltaics and energy-efficient ambient 
lighting mm- Thanks to the recent developments of adaptive wavefront shaping and phase 
recording techniques in optics, the spatial degree of freedom of the input light can be controlled 



at an unprecedented level of precision. These developments have enabled coherent control of 
light propagation in highly scattering media by manipulating the interference of multiply scat¬ 
tered waves 171644761 . One striking interference effect that has caught much attention is the exis¬ 
tence of highly transmitting channels, termed ’’open channels” in a diffusive system. These open 
channels, which enable an optimally prepared coherent input beam to transmit through a strong 
scattering medium with order unity efficiency, were predicted initially for electrons EM). 
Since it is much more difficult to control the input electron states than the input states of clas¬ 
sical waves, the wavefront shaping technique has been utilized in the past few years to increase 
the coupling of the incident light to the open channels of random media li6lf8l]9l 17744731 . The 
open channels greatly enhance light penetration into the scattering media, that will have a pro¬ 
found impact in a wide range of applications. 

In reality absorption exists in any material system, and could have a significant impact on 
light transport in both diffusion regime and localization regime (2444261 . On one hand the in¬ 
terference effects may be modified by absorption, on the other hand light absorption in strong 
scattering media can be drastically enhanced or suppressed by interference effects 12714311 . 
Thus the interplay between absorption and interference determines not only the amount of en¬ 
ergy being transmitted, but also the amount of energy being deposited in a random medium. 
For example, wavefront shaping has enabled focusing of laser light onto a localized absorber 
that is buried within a random medium to enhance the local absorption (32]. 

Absorption also has a strong effect on the transmission eigenchannels, especially the open 
channels. The transmission channels are eigenvectors of the matrix t ’ t, where t is the field 
transmission matrix of the system. The eigenvalues T are the transmittance of the corresponding 
eigenchannels. In the lossless diffusion regime, the density of the eigenvalues T has a bimodal 
distribution, with one peak at T ~ 0 that corresponds to closed channels, and a peak at T ~ 1 
that corresponds to open channels 17714231 . An open channel has a spatial profile extended 
throughout the entire random medium, with the intensity maximum near the center f33l . When 
strong absorption is introduced uniformly across the entire system, the diffusive transport of 
light in the maximum transmission channel turns into quasi-ballistic l26l . The straightening of 
optical paths through the random medium reduces the dwell time and minimize the attenuation 
by absorption. The statistical distribution of transmission eigenvalues are no longer bimodal, as 
the peak at T ~ 1 is diminished by strong absorption (24]. Experimentally absorbers are often 
distributed non-uniformly in random samples, and it is not clear how the open channels would 
respond to spatially inhomogeneous absorption. 

In this paper, we present a numerical study on the transmission eigenchannels in disordered 
waveguides with spatially localized absorbing regions. We calculate the statistical distribution 
of transmission eigenvalues and find it can be broader than that with uniform absorption. The 
longer tail of the distribution is attributed to the higher transmission eigenchannels that manage 
to circumvent the local absorbing regions to minimize loss. Compared to other transmission 
eigenchannels, the one with the largest eigenvalue is the most efficient in bypassing the ab¬ 
sorbing regions to transport the maximal amount of energy through the random system. Con¬ 
sequently, the ratio of the maximum transmittance to the average transmittance increases with 
absorption and eventually exceeds the ratio without absorption. 

2. Numerical model 

In our simulation, we consider a 2D disordered waveguide, shown schematically in Fig. 1(a). 
The dielectric cylinders with refractive index n = 2.0 and radius r c = 50 nm are randomly po¬ 
sitioned inside a waveguide with perfectly reflecting sidewalls. The dielectric cylinders occupy 
an area fraction of 0.04 corresponding to an average distance between cylinders of a = 0.44 11 m. 
The probe light enters the waveguide from the left open end and is scattered by the cylinders. 


The wavelength of input light A is set to 510 nm, to avoid the Mie resonances of individual 
dielectric cylinders. The light is transverse magnetic (TM) polarized, its electric field is parallel 
to the cylinder axis (z-axisj. The width of the waveguide is W = 10.5 ,um; the number of guided 
modes in the empty waveguide is IV = 2 W /A = 41. The length of the random array of cylinders 
is L =20.6 jl m. 



Fig. 1. Schematic of the 2D disordered waveguide used in our numerical simulation. Di¬ 
electric cylinders are placed randomly in a waveguide with perfect-reflecting sidewalls, (a) 

Qn: homogeneous distribution of absorbers across the entire random arrays of cylinders. 

(b) < 23 : absorbers are confined to the three circles with diameter D a . (c) Q \: absorbers are 
concentrated in a single circular region of diameter D u in the middle of the disordered 
waveguide. 

We use the recursive Green’s function method H291l341f35l to calculate the transmission ma¬ 
trix of the disordered waveguide |[26l . which gives the output field for any arbitrary input. 
Using the input and output fields as the boundary conditions, we further compute the field dis¬ 
tribution inside the disordered waveguide. The field intensity is averaged over the waveguide 
cross-section to give the evolution I(x) along the waveguide (in the x direction). The ensemble- 
averaged (!(x)) displays a linear decay, from which we extract the transport mean free path l, = 
1.65 fim f36| . The localization length is then estimated to be t, = (n/2)Nl t =107 /im. Since 
the propagation of light in the disordered waveguide can be described by diffusion. 

After characterizing the scattering properties, we introduce optical absorption to the disor¬ 
dered waveguide. Below we consider three cases. The first one, labeled Q u , has homogeneous 
absorption across the entire random structure [Fig.QJa)]. More specifically, a constant imag- 










inary refractive index, «,■ > 0, is introduced to both dielectric cylinders and background, to 
avoid any additional scattering caused by the spatial inhomogeneity of m. The ballistic absorp¬ 
tion length is l a = 1/(2 kni) = \/(pO a ), where k = 2n/X. is the wavevector, p is the density 
of absorbers, and O a is the absorption cross-section of each absorber. The diffusive absorp¬ 
tion length, given by Q, = yjl t l a /2 , determines the strength of absorption effects. In the weak 
absorption regime, q„ > L, the average length of diffusive paths inside the random medium 
l p = 2L 2 /l, is shorter than l a , thus most scattering paths are barely affected by absorption. In 
the strong absorption regime, q a < L, large attenuation of long scattering paths significantly 
modifies the transport through the system j26l . 

The next two cases have non-uniform absorption in the disordered waveguides, one is labeled 
Q\ where the absorption is confined to three isolated circles [Fig. 0b)], the other is Q\ which 
has a single absorbing region in the middle of the waveguide [Fig. Q](c)]. All the circular 
absorbing regions have the same diameter D a = 2.7/,, and the space in between the absorbing 
regions, as well as the distance from an absorbing region to the waveguide sidewall, is larger 
than the transport mean free path. When comparing the effects of absorption in the above three 
cases, the total number of absorbers is kept constant, so that only the spatial distribution of 
absorbers is different. Since n, is proportional to the density of absorbers p, the product iijS a 
is the same, where S a is the total area of absorbing regions. The smaller the S a , the larger is 
the value of n,-, and the stronger is the absorption within the absorbing region. Hence, l a °c S a 
and °c y/Sj,. Below, the absorption strength is given by aL/q a , where a = y/S a /(LW) is a 
scaling factor which is equal to 0.27 for Q\, 0.46 for Qy and 1.0 for Q u . 


3. Statistical distribution of transmission eigenvalues 

A singular value decomposition of the transmission matrix t gives 

t = UHV\ (1) 


where E is a diagonal matrix with non-negative real numbers, o n = \f%,, T„ is the transmittance 
of the n th transmission eigenchannel, i\ > T2 > T3... > T,y. U and V are N xN unitary matrix, 
V maps input channels of the empty waveguide to eigenchannels of the disordered waveguide, 
and U maps eigenchannels to output channels. The input singular vector that corresponds to the 
highest transmission eigenvalue T\ gives the maximum transmission eigenchannel, its elements 
represent the complex coefficients of the waveguide modes that combine to achieve the highest 
transmission through the random waveguide. 

In the absence of absorption, the density of the transmission eigenvalues of the disordered 
waveguide P( t) [Fig. [2a)] follows the bimodal distribution [dashed line in Fig. [2a)]: fl7]423ll 
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where f represents the average transmittance. As shown in Fig. 2(a), the distribution has two 
peaks, one at x ~ 1 and another at x ~ 0. The transmission eigenchannels with T ~ 1 are “open 
channels”, and the ones at T ~ 0 are “closed channels”. The diffusive transport is dominated by 
the open channels, and f is determined by the number of open channels I1171I18II . Figure [2b) 
plots the spatial distribution of field amplitude | E z (x,y) \ for the maximum transmission channel 
in one random realization of the disordered waveguide. The input light penetrates through the 
entire waveguide. The cross-section-averaged intensity I(x) = (1 /W) \E z (x,y)\ 2 dy, shown 

in Fig. 2c), is peaked at the center of the waveguide l33l . 

When absorption is introduced uniformly across the disordered waveguide, the open channels 
experience more attenuation than the closed channels since light in an open channel propagates 








Fig. 2. Statistical distribution of transmission eigenvalues P(z) for disordered waveguide 
with homogeneous and inhomogeneous absorption, (a) Without absorption, P( z) (filled 
circles) exhibits the bimodal distribution given by Eq. (2) (dashed line), (b) Spatial dis¬ 
tribution of the electric field amplitude |£j(jr,y)| of the highest transmission eigenchannel 
in one random realization of disordered waveguide, (c) Cross-section-averaged intensity 
along the x-direction, having the maximum at the center of the waveguide, which coincides 
with the absorbing region in Q\ (marked by the vertical dashed lines), (d) P(z) for weak ab¬ 
sorption aL/^a = 0.4 [a = 1 for Q u , a = 0.46 for Qt, and a = 0.27 for Q\\. The peak near 
T = 1 is diminished and shifted towards smaller z. The shift in the case of inhomogeneous 
absorption is slightly larger than that of homogeneous absorption, as the open channels 
experience more attenuation due to better spatial overlap with the localized absorbing re¬ 
gions. (e) P(z) for strong absorption aL/= 1.8. Q\ has the highest transmission among 
the three cases, (f) The normalized width W t of P(z) as a function of the absorption strength 
aL/^ a . (g, h, i) Spatial map of the normalized Poynting vector S' (x,y) for the maximum 
transmission channel in the disordered waveguide with homogeneous absorption Q u (g), 
three absorbing regions 03 (h) and one absorbing region Q\ (i) at aL/£, a = 1.8. Dashed 
circles in (h, i) mark the boundary of the absorbing regions. When absorption is strong and 
inhomogeneous, main energy flows bypass the absorbing regions. 




















deeper into the waveguide. Consequently, the peak of P( t) near T = 1 is diminished and shifted 
towards smaller T [dotted line in Fig.(2jd)]. The effect of absorption can be further enhanced by 
inhomogeneous absorption, e.g. in the case of Qi, all absorbers are concentrated in the center of 
the waveguide, which coincide with the intensity maximum of the highest transmission channel 
[dashed line in Fig. 2(c)]. The calculated P( t), shown in Fig. 2(d), has the T — 1 peak moved 
farther to smaller T, reflecting a faster decreasing transmission of the open channels in Q\. 

However, as we continue increasing absorption, the behavior changes completely. As seen 
in Fig. |2e), the peak of P(z) at large T vanishes in all three cases, and the bimodal distribu¬ 
tion is replaced by a monotonic decay of P( t) with T. Surprisingly, the disordered waveguide 
with a single absorbing region Q\ exhibits the longest tail at high transmission, followed by 
the waveguide with three separate absorbing regions {?3, while the waveguide with uniform 
absorption Q u has the shortest tail. This trend is just opposite to that with weak absorption. To 
quantify the change in the width of P( t), we plot the normalized width W r = 1 /(t 2 )/(t ) 2 — 1 
in Fig. |2f). While W T for Q„ decreases monotonically with increasing absorption, W z for Q 3 
and Qi first decreases and then increases, eventually exceeds the value without absorption. 

4. Maximum transmission channel 

To understand why inhomogeneous absorption can reach higher transmission than homoge¬ 
neous absorption, we examine the maximum transmission channel in the presence of strong 
absorption. To map the energy flow inside the disordered waveguide, we compute the Poynting 
vector S(x,y) = (1 /2)Re[£ (x,y) x H* (x,y)]. The net flow over a cross-section of the disordered 
waveguide is F(x) = Jq S x (x,y)dy, where S x (x,y ) is the projection of S(x,y) on the x-axis. 
While the net flux F(x) remains constant in the absence absorption, it decays in the presence 
of absorption. For a clear visualization of the energy flow deep inside the random structure, 
the Poynting vector S(x,y) is normalized by F(x), S'(x,y ) = S(x,y)/F(x). Figure EJg-i) plot 
the normalized Poynting vector S'(x,y ) for the maximum transmission eigenchannels in the 
disordered waveguide with homogeneous or inhomogeneous absorption. With homogeneous 
absorption Q u [Fig. |2g)], light propagates mostly in the forward direction in order to reduce 
the dwell time within the random waveguide to minimize loss. Light transport changes from 
diffusive to quasi-ballistic when absorption is strong (26]. In disordered waveguide with three 
separate absorbing regions Qj [Fig.[2] h)], the main energy flows are in between the three sepa¬ 
rate absorbing regions to avoid the absorption. When there is only one absorbing region in the 
middle of the waveguide Q\ [Figure^)], the incoming energy stream splits into two to circum¬ 
vent the absorbing region. Therefore, the high transmission channels have modified their energy 
flows to bypass the absorbing regions in the waveguide so as to achieve higher transmission in 
the case of inhomogeneous absorption than that of homogeneous absorption. 

Next, we track the evolution of the highest transmission channel with a gradual increase 
of absorption. Figure [3 a) plots the change in the maximum transmission eigenvalue Tf as a 
function of the absorption strength ttL/i^ a for the disordered waveguide with homogeneous or 
inhomogeneous absorption. When absorption is weak aL/£, a < 1, Ti for Q \ decreases faster 
than Qt, and Q„. However, as absorption becomes strong aL/t, a > 1, the drop of Ti slows 
down and its value for Q\ is higher than the other two cases. This transition occurs when the 
maximum transmission channel modifies its energy flow to bypass the absorbing region. For 
comparison, we plot the average of transmission eigenvalues f as a function of absorption 
strength in Fig. [3jb). Its decay is qualitatively similar to that of T \, because f has the largest 
contribution from Z \. Their ratio T] /f displays the subtle difference in the change of their decay 
rates by absorption. In the weak absorption regime [Fig. [3c)], the ensemble-averaged ratio 
(Ti /f) decreases monotonically for all three cases. The reduction is the fastest for Q \ when the 
number of total absorbers is the same. This can be understood from the results in the previous 



<Ti/T ) 



Fig. 3. Comparison of the maximum transmission eigenvalue tj and the average of all trans¬ 
mission eigenvalues f in the presence of inhomogeneous absorption to that of homogeneous 
absorption, (a) T\ decreases with absorption in all three cases. The drop is the fastest for 
Q\ when absorption is weak [aL/^ a < 1] but switches to the slowest when absorption is 
strong [ocL/£ a > 1], (b) Similar trends are observed for the change of f with absorption, 
(c) The ensemble-averaged ratio (tj /f) shows the fastest reduction for <2i compared to Q u 
and <23 in the weak absorption regime, (d) (ti/t) starts to increase in the strong absorption 
regime, and eventually exceeds the ratio without absorption in Q\ and <23• 













section. The spatial distribution of the maximum transmission channel is barely modified by 
weak absorption, and it experiences the most attenuation when all absorbers are concentrated 
in the spatial location where its field intensity is maximal ( Q \). The lower transmission channels 
do not penetrate as deep into the random waveguide, thus their intensity maxima are closer to 
the input end of the wavewguide. The less spatial overlap with the absorbing region leads to 
lower attenuation, thus the average of transmission eigenvalues x decreases more slowly than 
H- 

With a further increase of absorption, the ratio X\fx starts increasing, especially in the case of 
inhomogeneous absorption (Q 3 , Q\), its value eventually exceed that without absorption. This is 
because the maximum transmission channel bypasses the absorbing regions to reduce the loss, 
and its transmittance decreases less than other eigenchannels. Thus the reduction of X\ becomes 
smaller than that of f, leading to an increase of their ratio Ti/x with absorption. The ratio X\ /i 
surpasses its value without absorption at D a /C, a ' V! « 0.6 for (F, and D a /^a ' 1 « 0.8 for Q\, where 
4 !and fjj 1 j are the diffusive absorption lengths within the absorbing regions in (9 3 and Q \ 
respectively. Hence, the highest transmission eigenchannel starts to circumvent the absorbing 
regions when the diffusive absorption length becomes comparable to the size of the absorbing 
region. In the case of homogeneous absorption, X\/x levels off and then increases slightly when 
aL/^ a > 1 [Fig-Eld)]- Since the absorbers are everywhere in the random waveguide, light in 
the maximum transmission channel cannot bypass absorbers when transmitting through the 
waveguide. The only way to minimize loss is to shorten the dwell time inside the waveguide 
by taking less winding paths when absorption is strong. The other transmission channels do not 
change as much as the maximum transmission channel, so they may experience slightly higher 
loss. 

5. Other transmission eigenchannels 

In this section, we investigate the changes in other transmission eigenchannels due to inho¬ 
mogeneous absorption. Figure |4j a) shows the ratio t„/t° for Q \ , where z n is the transmission 
eigenvalue for the n' h channel with absorption and is the eigenvalue without absorption. 
When absorption is weak ' < 1, Tf decreases faster than all others. We also calculate the 

absorption in each channel, A n = 1 — x n — R„, where R„ is the reflectance of the n ,h transmission 
eigenchannel. The absorption in the maximum transmission channel A\ is the largest [Fig.[4]b)j, 
which explains why X\ reduces more than other channels when absorption is weak. However, 
at higher absorption level D a /^a 1 > 1 , the transmission reduction of other eigenchannels such 
as T 3 exceeds that of X\ . Correspondingly, the absorption experienced by the third eigenchannel 
A 3 becomes larger than the first eigenchannel A 1 [Fig. |4jb)]. 

The above result can be explained by the modification of these transmission channels by in¬ 
homogeneous absorption. In a lossless disordered waveguide, the field intensity maximum of 
the first transmission eigenchannel is located at the center of the waveguide, whereas the inten¬ 
sity maximum of the third eigenchannel is shifted towards to the input end of the waveguide 
[Fig. He, d)]. When absorption is introduced only to the central region of the waveguide, the 
field intensity of the first channel is peaked in the middle of the absorbing region while the field 
intensity of the third channel is peaked at the edge [black dashed lines in Fig. [4jc, d)]. With 
weak absorption 11 > D a , the spatial distribution of field intensity in both channels are barely 
changed, thus the first eigenchannel is attenuated more than the third eigenchannel [brown dot¬ 
ted lines in [4jc . d)]. However, at strong absorption where the diffusive absorption length within 
the absorbing region becomes shorter than the size of the absorbing region < D a , light 
transmission through the central absorbing region becomes very low and the scattering paths in 
the first eigenchannel avoid that region to achieve high transmission as seen in Fig. [2ji)- The 
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Fig. 4. Modification of transmission eigenchannels in disordered waveguide with a single 
absorbing region in the middle Q\. (a) Ensemble-averaged ratio (t„/t®) for individual 
eigenchannels. The reduction of Tj is the largest when absorption is weak D a /^i ^ = 0.32 
but becomes less than other eigenchannels such as T 3 when absorption is strong D a /^ 1 ' 1 = 
1.45. (b) The absorption of individual eigenchannels A„. The first eigenchannel is the one 
that experiences the most absorption at weak absorption but it is replaced by the third 
eigenchannel at higher absorption level, (c, d) Ensemble average of cross-section-averaged 
electric field intensity (7(x)} for the first and third eigenchannels. Both are normalized to 
their maximal values when there is no absorption. With weak absorption D a = 0.32 
(brown dotted lines), the field intensity of the first eigenchannel decreases more at the 
waveguide center than the third eigenchannel. The vertical dashed lines mark the boundary 
of the absorbing region. When absorption is strong D a f tj u 1 ^ = 1.45, the reduction of field 
intensity behind the absorbing region becomes less for the first eigenchannel than the third 
eigenchannel as the former is modified to circumvent the absorbing region. 







circumvention of light paths around the absorbing region helps to reduce the absorption of light 
and therefore the field intensity in the section of the waveguide behind the absorbing region is 
higher for the first eigenchannel than the third eigenchannel [green solid lines in Fig.|4jc, d)]. 

Similar modifications of the transmission eigenchannels are observed in the disordered 
waveguide with three separate absorbing regions Qj. Since it is more difficult to bypass three 
absorbing regions than a single one, the decrease of X\ is higher for (F, than for Q\ [Fig. [3 a)], 
and also Qt, has a narrower distribution of transmission eigenvalues than Q \ [Fig. |2je,f)] . Note 
that in both Q\ and Q 3 , the space in between the absorbing regions or between the absorbing 
region and the waveguide wall is larger than the transport mean free path, so that the multiple 
scattering of light in the non-absorbing regions and the interference of the scattered light enable 
an efficient steering of energy flow away from the absorbing regions. 

6. Conclusion 

We have performed a detailed numerical study to understand how spatially non-uniform ab¬ 
sorption modifies the transmission eigenchannels in a 2D disordered waveguide. In the weak 
absorption regime, the spatial profile of a transmission channel remains very similar to that 
without absorption, and the effect of inhomogeneous absorption can be stronger or weaker than 
homogeneous absorption depending on the spatial overlap of the localized absorbing region 
with the field intensity maximum of the channel. In the strong absorption regime, the high 
transmission channels redirect the energy flows to circumvent the absorbing region to mini¬ 
mize loss. Thus the attenuation of high transmission channels by inhomogeneous absorption 
is lower than that by homogeneous absorption, making the statistical distribution of transmis¬ 
sion eigenvalues in the former broader than that in the latter. Compared to other transmission 
eigenchannels, the one with the largest eigenvalue is the most efficient in bypassing the absorb¬ 
ing regions to transport the maximal amount of energy through the random medium. Hence, 
the ratio of the maximum transmittance to the average transmittance, as well as the normalized 
width of the eigenvalue distribution, increases with absorption and eventually exceeds the value 
without absorption. Our numerical study provides a physical understanding of the effects of in¬ 
homogeneous absorption on transmission eigenchannels in diffusive media. The finding that 
inhomogeneous absorption may have a weaker impact on open channels than homogeneous 
absorption is promising for practical applications. 
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